Magnetotransport properties of a magnetically modulated two-dimensional electron 

gas with the spin-orbit interaction 
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We study the electrical transport properties of a two-dimensional electron gas with the Rashba 
spin-orbit interaction in presence of a constant perpendicular magnetic field (Bqz) which is weakly 
modulated by Bi = B 1 cos(qx)z, where Bi <C Bo an d q = 2-7r/a with a is the modulation period. 
We obtain the analytical expressions of the diffusive conductivities for spin-up and spin-down elec- 
trons. The conductivities for spin-up and spin-down electrons oscillate with different frequencies 
and produce beating patterns in the amplitude of the Weiss and Shubnikov-de Haas oscillations. 
We show that the Rashba strength can be determined by analyzing the beating pattern in the 
Weiss oscillation. We find a simple equation which determines the Rashba spin-orbit interaction 
strength if the number of Weiss oscillations between any two successive nodes is known from the 
experiment. We compare our results with the electrically modulated 2DEG with the Rashba inter- 
action. For completeness, we also study the beating pattern formation in the collisional and the 
Hall conductivities. 

PACS numbers: 71.70.Ej,73.43.Qt,85.75.-d 



I. INTRODUCTION 

The magnetotransport properties of a two-dimensional 
electron gas (2DEG) in presence of a weakly modulated 
one-dimensional (ID) periodic electric potential has been 
studied in great details experimentally and theoretically 
for a long timei~— . In the absence of the modulation 
and the spin-orbit interaction (SOI), the magnetocon- 
ductivity oscillations due to the charged impurities are 
commonly known as the Shubnikov-dc Hass (SdH) oscil- 
lations. In the presence of the modulation and at low 
magnetic fields, it has been observed that the magne- 
toresistivity tensor oscillates with the inverse of the mag- 
netic field. These oscillations are completely different in 
periodicity and temperature dependence from the SdH 
oscillations observed at higher magnetic field. These os- 
cillations are commonly known as the Weiss oscillations. 
This is due to the effect of the commcnsurability between 
the two length scales in the system: the cyclotron diam- 
eter at the Fermi energy and the modulation period a. 
The similarity and differences between the Weiss and the 
SdH oscillations are as follows: i) both the oscillations 
are periodic in 1/B ii) the period of the Weiss oscilla- 
tion varies with the electron density (n e ) as y/rQ, whereas 
that of the SdH ones as n e ; iii) the amplitude of the Weiss 
oscillations less dependent on temperature than that of 
the SdH oscillations; iv) the Weiss oscillations are visible 
at weak magnetic field (Bq < 0.6 T) and the SdH oscil- 
lations modulated by the Weiss oscillations are visible at 
higher fields. 

In conventional 2DEG systems, magnetotransport 
properties in the presence of ID electric and magnetic 
modulations are continuing to be an active research field. 
The magnetotransport properties of a 2DEG in presence 
of a weakly modulated magnetic field has been studied 
theoretically 6 - - — . It was theoretically observed that the 
magnetoconductivity oscillates with inverse of the mag- 



netic field and amplitude of the oscillation is much higher 
compare to the case of electrical modulation. The mag- 
netothermodynamical properties of electrically or mag- 
netically modulated 2DEG has also been studied where 
the Weiss-type oscillation is also showr* 5 ^ 13 ' 14 . Later, 
the magnetoresistance oscillation in a magnetically mod- 
ulated 2DEG system was observed experimentally-i 5 - - — . 
In these experiments, a ID spatially varying magnetic 
modulation was achieved by placing micro-patterned fer- 
romagnet or superconductor on the surface of a 2DEG 
system. 

An internally generated crystal field induces the SOI 
which is known as the Rashba SOI. The Rashba inter- 
action strength can also be controlled by a strong ex- 
ternal electric field acting normal to the 2DEG plane. 
Using the Rashba interaction it was possible to ex- 
plain many experimentally observed features like the 
combined resonance a 18 ' 19 and the beating patterns in 
the SdH oscillations 20 . The Rashba interaction is re- 
sponsible for many other novel effects like the spin- 
FET— , spin-galvanic effect^ and the spin Hall effec1s22i2i. 
The Rashba SOI in the 2DEG influences various prop- 
erties such as transpor t 25 ' 26 , magnetotransport^! - — , 
magnetization^! etc. In Refsi 3 ^ - — , the effect of the 
Rashba interaction on magnetotransport properties of 
electrically modulated 2DEG have been studied theoret- 
ically. 

Generally, beating pattern analysis in the SdH 
oscillation^ - and weak anti-localization method^ 5 , are 
used to extract the strength of the Rashba SOI in the 
2DEG system. However, in the presence of Bq ^ 0, Zee- 
man splitting is also accompanied with the Rashba spin 
splitting. To avoid this problem optical measurement^ 6 - 
has also been proposed. The Weiss oscillation is due to 
the modulation induced Landau levels broadening which 
oscillates with inverse magnetic field. Moreover, this os- 
cillation appears at low magnetic field where the Rashba 
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SOI dominates over the Zeeman splitting. The Weiss os- 
cillation is less influenced by the Zeeman term specially 
in the periodicity. To determine the Rashba strength, 
beating pattern analysis in the Weiss oscillation is more 
reliable than the SdH oscillation. 

In this work, our primary goal is to study analytically 
the Weiss oscillations of the 2DEG with the Rashba in- 
teraction and to determine the strength of the Rashba 
SOI by analyzing the beating pattern of the Weiss oscil- 
lations. For completeness, we also study other transport 
coefficients such as the collisional and the Hall conduc- 
tivities of the spin-up and spin-down electrons in details. 

This paper is organized as follows. In section II, 
we summarize the energy eigenvalues, the correspond- 
ing eigenfunction and density of states of the 2DEG with 
the Rashba interaction in presence of a uniform magnetic 
field. Also, the first-order energy corrections due to the 
magnetic perturbations and group velocities for spin-up 
and spin-down electrons are evaluated. We calculate and 
discuss the diffusive, collisional and the Hall conductivi- 
ties by using the semi-classical Kubo formula in section 
III. The summary of our work is presented in section IV. 



II. ENERGY EIGENVALUE, EIGENFUNCTION 
AND DENSITY OF STATES IN PRESENCE OF 
THE RASHBA SOI 

We consider a 2DEG in the x-y plane subjected to a 
magnetic field B = [Bo + B\ cos(qa;)]z, where q — 2ir/a 
and a is the modulation period. Also, we consider 
strength of the magnetic modulation is very weak i.e., 
B\ <C Bo- In the Landau gauge the corresponding vector 
potential A is A = [Box + (Bi/q) s\\\(qx)]y. 

The Hamiltonian of an electron with charge — e in pres- 
ence of the perpendicular magnetic field B is 

H = (P ^f )2 l+|[<rx(p + eA)] z + | Mb B ■ a, (1) 

where p is the 2D momentum operator, m* is the effective 
mass of the electron, g is the Lande g- factor, [i B is the 
Bohr magneton, 1 is the identity matrix, <r = (er x , a y ,a z ) 
are the Pauli spin matrices, and a is the strength of the 
Rashba interaction. 

Expanding the above Hamiltonian and re-writing as a 
sum of the various Hamiltonians: H = H + Hi + H 2 + 
H 3 + H4, where 

1 a 
H = \Px + (Py + eB x) 2 ]t + - [<r x (p y + eB x) 

~ °yPx] + ^9V B B cr z , (2) 



TT 

Hi 


= ~^~(Py + eBix) sin(qx)±, 


(A) 


Hi 


= V n 7r-<Jx sin(ga;), 

Z7T 


(4) 


H 3 


= V B g* cos(qx)a z , 


(5) 


Hi 


= ^K {l ~ C0<2qx)}1 - 


(6) 



Here, V B = hui = heBi/m* is the strength of the 
effective magnetic potential determined by the ampli- 
tude B\ of the magnetic modulation and uj\ <C ujq with 
ljo = eBo/m* is the cyclotron frequency due to the con- 
stant magnetic field Bo- The dimensionless parameter 
p a = ak a with k a = am* /h 2 and g* — gm* / (4mo). 
Also, e a = h 2 /(m*a 2 ) is a characteristic energy scale in- 
troduced by the modulation period a. Here, the term Hq 
is the Hamiltonian for the electron with the Rashba in- 
teraction in presence of the constant magnetic field Bo z 
including the Zeeman energy. The terms H\ and H4 are 
due to effect of the magnetic modulation on the kinetic 
energy term. The term Hi is the interaction between 
the spin-orbit and the magnetic modulation whereas the 
term H3 is the interaction between the Zeeman energy 
and the magnetic modulation. 

The Hamiltonian Ho can be solved analytically and 
can be treated as a unperturbed Hamiltonian. The other 
four Hamiltonians can be treated as a small perturba- 
tions since V B <C ftioo- The eigenfunctions of the unper- 
turbed Hamiltonian (Ho) can be used to find the energy 
correction due to the small perturbations Hi , Hi , H 3 and 
H A . 

Here, we shall briefly summarize the results of the 
Refill where the analytical solutions of the Hamiltonian 
Hq have been derived. Using the Landau wave func- 
tions without the Rashba interaction as the basis, one 
can write new wave function as 

ik v y °° / r< + \ 

* fc9 (r) = K{x + xo) ~» . (7) 

V L V n=0 V L/ « / 

Here, <p n (x) = (1/ ^ ^2 n n\lo)e~ x2 / 2l °H n (x/lo) is the 
normalized harmonic oscillator wave function with n is 
the Landau level index, L y is the width of the sample in 
the y direction, Iq = ^Jh/(eBo) is the magnetic length, 
and the cyclotron orbit is centered at — xq with xq = kyly. 

Using these wave functions and Eq.(l), the eigenvalue 
problem Ho^ — E^ff leads to an infinite number of equa- 
tions that can be solved exactly after decomposing it into 
independent systems of one or two equations 2 —. The re- 
sulting eigenstates are labeled by a new quantum num- 
ber s instead of n. For s = 0, there is only one level, the 
same as the lowest Landau level without SOI, with en- 
ergy Eq = Eo = (hwo — gfi B Bo)/2 and the corresponding 
wave function is 

*+ Jks (r) = -^^ (x + x )^J. (8) 
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For s = 1,2,3.... there are two branches of energy levels, 
denoted by + and — with energies 



ET = shu}„± x /E 2 + sEJvjj q , (9) 



where E a — 2m*a 2 /h 2 . These energy levels are degener- 
ate in k y . Therefore, the group velocities, v x and v y , are 
zero. 

The corresponding wave function for + branch is 



y/LyA a 

and for — branch is 

p ikyy 

$- ( r ) = 

^ ' ^Ly-A S 



D s 4> s ^i(x + x ) 

(j> 8 {x + X ) 



<ps-l(x + Xq) 
-D s 4> s (x + Xq) 



(10) 



(11) 



where A s = 1 + D 2 and D s = y/sE a haj /[Eo + 
V 'El + sEJiwo]. 

The Landau level quantum numbers s at the Fermi 
energy are determined from the equation E F = sHuiq ± 
\J E 2 + sE a frwo- These quantum numbers are given 
by s^ 1 ~ k 2 l 2 /2^fk a k F l 2 , where k F — y / 2nn e is the Fermi 
wave vector and n e is the electron density. 

The approximate density of states (DOS) of the 2DEG 
with the Rashba SOI in the presence of constant perpen- 
dicular magnetic field is given 30 by 



D ± (E) 



2ttH 2 



2 exp | 



( hujQ 



r 2tt / E a 



E„E 



,12) 



where To is the Landau level broadening. The detail 
derivation of the DOS is given in the Appendix A. We 
shall use the above DOS to obtain the analytic expres- 
sions of the magnetotransport coefficients in the presence 
of the modulations. 

Using the perturbation theory, we calculate the first- 
order energy correction for + branch as well as — branch. 
It is to be noted that the correction due to the Hamilto- 
nian is of the order of lu 2 which is much smaller than 
ujq. Therefore, we neglect the contribution of H<± to the 
total energy correction. 

The total energy for + branch and — branch is given 

by 



E tk y = shu)a±^El + sE a hu + cQs(qxo)F±(ull3) 
where 

FHu) = ^[e-^g*{±D 2 s L s _ 1/2 ^ /2 (u) 
T L s _ 1/2 ± 1/2 (u)\ + D 2 G a _ l/2 ^i/2 (u) 
+ G s _ 1/2±1/2 (u) ± -|=e-»/ 2 Li_ 1 («)].(14) 
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FIG. 1: Plots of the bandwidth A ± /Vb at the Fermi energy 
versus dimensionless inverse magnetic field B a /Bo for two dif- 
ferent values of a. Here, solid and dashed lines correspond to 
the upper branch and lower branch, respectively. 



Here, u = q 2 Iq/2, s — 0,1,2... for + branch and s = 
1,2,3... for - branch. Also, G s (u) = e - u / 2 [L s (u)/2 + 
Ll_i(v)] and Z = 2k a DJo. The upper and lower signs 
correspond to the + and — branches, respectively. The 
width of the broadened levels of the two branches is |A^| 
and it can be written as A^ = 2F^(u). 

For various plots, we use the following parameters: a 
in units of ao = 1.0 x 10 -11 eV-m, modulation strength 
V B = 0.05 meV for Bi = 0.02 T, electron density n e = 
10 16 /m 2 , electron effective mass m* = 0.05mo with mo is 
the free electron mass, mobility fi = 100 m 2 /V-s, g = 2, 
modulation period a — 800 A and temperature T = 1.5 
K. For these parameters, p a — 1.05 for a — 2ao, p F = 
ak F = 20, E F = 47.5 meV and e a = 0.236 meV. 

In Fig. 1, we plot the dimensionless bandwidth A ± /Vb 
at the Fermi energy as a function of the dimensionless 
inverse magnetic field A = B a /Bo for two different values 
of a. Here, B a = h/(ea 2 ) = 0.102 T is the characteristic 
magnetic field introduced by the modulation period a. 
We plot A^ 1 instead of | A ± | to show the oscillations more 
clearly. 

The diagonal matrix elements of the velocity operator 
do not vanish due to the k y dependence of the energy lev- 
els. It can be calculated by using 
and their values are given by 



m)dEf /dky, 



41 , 



sm{qx )F s (u). 



(15) 



These non-zero values give rise to the finite diffusive con- 
ductivity whereas in the absence of the magnetic mod- 
ulation the diffusive conductivity vanishes whether the 
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Rashba interaction is present or not. These velocities for 
spin up and spin down branches do oscillate with slightly 
different frequencies as a function of the magnetic field 

So- 



ul MAGNETOTRANSPORT COEFFICIENTS 

In the presence of the weak modulation, there are two 
contributions to the transport properties: diffusive and 
collisional contributions. The former is due to finite drift 
velocity gained by the electrons in the presence of the 
modulation. The latter is the contribution from the hop- 
ping of the localized states due to scattering by impuri- 
ties. The diffusive conductivity decreases with increasing 
impurity scattering whereas the collisional conductivity 
increases with the increase of the impurity scattering. 
The collisional scattering is dominant at low tempera- 
ture. 

We follow the formulation of Refill to calculate the 
magnetotransport coefficients in the presence of the mod- 
ulation. In the linear response regime and weak scat- 
tering potentials, the conductivity tensor, in the one- 
electron approximation, is a sum of a diagonal and non- 
diagonal terms i.e. ct m „ = ct^„ + cr"jj, where — x,y. 
Generally, the diagonal conductivity contains both dif- 



fusion and collisional contributions: ct^„ = ct~; , 
whereas the Hall conductance is obtained from the non- 
diagonal conductivity cr^. In our problem, a xx = CT™ 1 

because a xx = and a yy = ct™ 1 + a yy . The Weiss os- 
cillation is observed in the diagonal component of the 
magnetoresistance p xx which is inverse of the conductiv- 
ity tensor: p x 



xy 



the modulation in the transport coefficients, ct™ 1 and a, 
is very weak. Therefore, the change in the conductivities 
due to the modulation [Act^j, = ct^^Vb) — ct^^O)] as a 
function of the magnetic field will be shown. 



A. Diffusive conductivity 

The standard semi-classical expression for the diffusive 
conductivity is given as22, 



dif 



So 



(16) 



where So is the area of the system, r is the electron relax- 
ation time, f{E) is the Fermi-Dirac distribution function, 
and j3 = (l/k B T) is the inverse of the thermal energy. 
Also, £ is the set of the quantum numbers: £ : {s, ct, k y } 
with a = ± and i& = ^ s the diagonal matrix 

elements of the velocity operator v y . 

The non-zero expectation value of the y-component of 
the velocity operator is reflected in the finite value diffu- 
sive conductivity whether the SOI is considered or not. 
In the presence of the Rashba interaction, the two energy 
branches, + branch and — branch, contribute to the total 



conductivity. The diffusive conductivity for ± branches 
is given as 



CT 



dif(±) 



/3e 2 



So 



E 



dk y f(Ef)[l-f{Ef)}{vi 



(17) 

In the limit of weak magnetic field Bo, many Landau 
levels are filled (s ^> 1), we use the following approxi- 
mations: e~ u / 2 L s {u) ~ (7T 2 sw) -1 / 4 cos(2- s /iu — tt/4) and 
E s -> Jo°° dE/{hwo). With the use of the above approx- 
imations and following Ref.^, Eq. (fTT|) reduced to an 
analytical form which is purely Weiss contribution, as 



Weiss 



-A B C±\ 



x sin (2nQ. ± X - 



(18) 



where il ± = 2(p F T P a ) are the frequencies of the con- 
ductivity oscillations of spin-up and spin-down electrons, 
A B = V 2 r/(he a ), T a — EF/(2ir 2 k B p F \) is the character- 
istic temperature, and $q = 5^ + Air/{p F T-Pa) are the 
phase factors. Other parameters are H(x) — x/ smh(x), 



4(p F TP J 



16tt 4 



k 



2tt 



1 



4{p F TP a ) 16^ 4 



■ sin 



( ** ). 

\P F T P r ' ' 



cr yy I \(T xx Uyy - cr xy a yx ). The effect of and 



± = (p,TPa) Bin /_2 !L _ 

4tt 2 V F TP„ 



with = ^{Cf) 2 + {C±) 2 . Also, 

tan -1 (Cf /Cf) + it. The amplitudes, frequencies and 
phase factors are different for different branches. Also, 
the amplitudes, frequencies and phase factors depend on 
n e , a and a. The diffusive conductivities for spin- up and 
spin-down electrons are oscillating with different frequen- 
cies in units of A. The total Weiss contribution to the 
diffusive conductivity is given by 



dif 

Weiss 



c 



H(^r){w+ sin (27rn+A-$+) 



W "sin (27rft _ A 



.">} 



(19) 



where C — Cq + Cq . It shows that the total diffusive 
conductivity exhibits beating patterns in the amplitude 
of the Weiss oscillations. In Fig. 2, we compare the ana- 
lytical result of the Weiss contribution (fT9|) to the total 
diffusive conductivity with the exact numerical results 
obtained from Eq. p7|) . The analytical result are in 
excellent agreement with the numerical result except at 
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higher magnetic field (A < 0.15) where the SdH oscilla- 
tion dominates over the magnetic Weiss oscillation. 

To obtain the analytical expressions of the locations of 
the beat nodes and the number of oscillations between 
any two successive nodes, we simplify Eq. (| 19[) by con- 
sidering p F ^> p a and l/(p F T Pa) — 0. Equation (fT9|) 
reduces to the following form: 



"Weiss 



" h 2tt 2 AbX 
x cos (2irVld\ 



1 



#(jr) sin (27rfi a A + tt) 



(20) 



where tt a = (fi + +Q-)/2 and n d = - fi_)/2. 



At the node positions Bq = Bj, we have the following 



condition: cos(27rAOd)|B =Bj = 0; which gives us 



8 Pa B a 

(2i + i) : 



(21) 



where j = 0,1,2,3,4.. are the j-th beat node and 
the corresponding magnetic field is Bj. By using the 
above equation, the beating nodes appear at B a /Bj = 
0.11,0.35,0.59,0.83,1.07... which are in excellent agree- 
ment with the exact numerical results. The above Eq. 
indicates that the magnetic field corresponding to the 
j = node is the lower limit of the Weiss oscillation 
below which the SdH oscillation starts to dominate. In 
practical, the numbering of the beat nodes is difficult. 
To remove this problem, the above equation can be re- 
written for any two successive beat nodes as 



1 

— =4B 



1 



B 



1 



(22) 



So, the strength of the Rashba SOI can be determined 
from the above equation by knowing the locations of the 
two successive beating nodes. The number of oscillation 
between any two successive beat nodes can be obtained 
from the first sine term of Eq. (|20|) . which is 



N osc = 2p F B a 



1 



B, 



1 



Using Eq. 



y+1 ' 
in the above equation, we get 



Pf 

2p a 



k F 
2k 



(23) 



(24) 



From the above equation, the following important con- 
clusions can be drawn, i) The above equation can be re- 
written as a = h 2 k F j (2m* N osc ) . Therefore, the Rashba 
SOI strength can be easily calculated by just counting the 
number of oscillation between any two successive nodes. 

ii) The number of oscillation between any two successive 
nodes is constant for given values of n e and a, whereas 
it depends on the magnetic field in the SdH oscillation 30 . 

iii) While the frequency of the Weiss oscillation depends 
on the modulation period but the number of oscillation 
between any two successive nodes does not depend on it. 
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FIG. 2: Plots of the exact (dashed) and asymptotic (solid) 
expressions of the diffusive conductivity versus dimensionless 
inverse magnetic field B a /Bo for electric and magnetic mod- 
ulations. Here strength of the Rashba SOI a = 2ao. 



The analytical expression of the diffusive conductivity 
given in Eq. (|19j) is not able to explain the origin of the 
superposition of the SdH oscillation on the Weiss oscilla- 
tion at higher magnetic field. This can be explained by 
using the asymptotic expression of the DOS given in Eq. 
(TT3J) and replacing the summation^ over discrete states 
by integration as J2 S —> 2nl 2 / °° D ± (E)dE, then we get 
dif(±) 



Vyy = Cweiss + °sdH> where ^Woiss S iven in Eq- CL 
and the Weiss oscillation modulated by the SdH oscilla- 
tion (og^jj) is given by 



SdH 



i^-{-«(S)K 



sin (2tt^ ± A 



x cos(2tt/ ± /B ), 



(25) 



where W ± = \J [Cf) 2 + (2C 2 ± ) 2 , $ ± = 
with S ± = tan" 1 (2Cf/Cf) + n. 



P. 



8 ± + 4tt/(p f t 
Also, /± = 



(m*/(fie)) [E F + E a /2 =p \/E 2 + E a E F \ are the frequen- 
cies of the SdH oscillations for spin-up and spin-down 
electrons in absence of the modulation^ and the char- 
acteristic temperature for the SdH oscillation is T c — 
hw /(2ir 2 k B ). For a = 2a n and T = 1.5 K, the 



ratio er 



dif 
SdH 



dif 

Weiss I 



is 0.3, 0.025, 0.001 for B a /B = 



0.3,0.5,0.8, respectively. Here, a 



dif 
SdH 



'SdH 



'SdH' 
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Note that the Weiss oscillation frequencies (Q^) for spin- 
up and spin-down electrons are different from the SdH 
oscillation frequencies / . 

To derive Eqs. (JT5J and we have used the DOS for 
unmodulated 2DEG. In the presence of the modulation, 
the DOS can be written as2£ a sum of the unmodulated 
(D ± {E)) and modulated (D±(E)) part as Df otBl (E) = 
D ± (E) + D±(E). Here, D ± {E) is given in Eq. {I5J) 
and (E) is of the order of Vg , which is very small in 
comparison to D ± (E). Moreover, the correction to the 
Weiss and the SdH oscillations due to the D m (E) will be 
of the order of Vg. The effect of the DOS correction due 
to the modulation on the Weiss and the SdH oscillations 
is really small and we have neglected it. 

Substituting a — in Eq. ([T5]) . we get the oscilla- 
tion period P = l/f2 = l/(2p F ) which is same as ob- 
tained in Ref£. We compare this theoretical result with 
the available experimental result^ for a = 0. In this 
experiment^, magnetic modulation period a = 1 /im and 
the electron density n e — 2.2 x 10 15 /m 2 . Using these pa- 
rameters, we obtain P = 0.00425 which is very close to 
the experimental result P ~ 0.00465. 

Comparison with the electrical modulation 
case: Here, we would like to compare the above men- 
tioned results for the magnetically modulated case with 
the electrically modulated system. The electrical modu- 
lation potential is described by H' = Ve cos^a;), where 
Ve is the amplitude of the electric modulation poten- 
tial. We have used Ve — 0.05 meV in our numerical 
calculation. In the case of electrically modulated system, 
numerical results have been discussed in Ref.—. Here, 
we provide an analytical expression of the diffusive con- 
ductivity for the electric modulation case. We obtain the 
Weiss and the SdH contributions to the diffusive conduc- 
tivity, which are given below: 



°Wciss 



A A 



~F P a 



■cos 2 {tt/0 f TP J} 



h (p. 

sin|27rrj ± A- 2vr/(p F =Fpj} 



l + H 



(26) 



Therefore, Eqs. ([22]). ((23)) and ([24]) for the magnetic mod- 
ulation case remain also valid in the electric modulation 
case. 

However, there are few important differences between 
the diffusive conductivities in the electric and magnetic 
modulation cases. There is a definite phase difference 
between the diffusive conductivities in electric and mag- 
netic modulation cases. The amplitude of the diffusive 
conductivity in the presence of the magnetic modulation 
is found to be nearly [p F /(27r)] 2 times higher in compar- 
ison to the electrical modulation case, which is the same 
as in absence of the Rashba SOI 6 . The role played by the 
Rashba SOI in the case of magnetic modulation is that, 
there is a difference between the amplitudes of conduc- 
tivity for the spin-up and spin-down electrons by a factor 
p Q /(27r 2 ). In the electric modulation case, amplitudes of 
the conductivities of spin up and spin down branches are 
nearly the same. 



B. Collisional Conductivity 

In electron systems the charge impurities play an im- 
portant role in magnetotransport properties. Collisional 
conductivity arises because of the migration of the cy- 
clotron orbit due to scattering from charge impurities. At 
low temperature, we can assume that electrons are elas- 
tically scattered by the charged impurities distributed 
uniformly. The standard expression for collisional con- 
ductivity is given by 3 -^ 



0$) 



(28) 



Here, = /j/ for elastic scattering, W^g is the transi- 
tion probability between one-electron states |£) and 
Also, = (£|r^|£) is the expectation value of the /z com- 
ponent of the position operator for the electron in state 
The scattering rate is given by 



and 



'SdH 



e 2 AX 

T7 — t — ^ cos 2 {tt I (p F TP a )} 
h (Pf T P a ) 

-{-»(£)>© 

1 + sin{2 7 rrj ± A - 2tt/( Pf t Pa)} 



x cos (27r/ ± /B ), 



(27) 



where A E = V F T/(hs a ). The lower panel of Fig. 2 shows 
the comparison between the analytical expression given 
in Eq. (|26l) and the numerical results reproduced from 
Eq. (20) of Ref. 32 . Our analytical expression matches 
very well with the numerical result. The Weiss and the 
SdH oscillation frequencies for spin-up and spin-down 
electrons are the same for both type of the modulations. 



q„-(r-R) 



\0\ 2 S(Es-E e ), (29) 



where q = q 0x x + q 0y y is a 2D wave- vector and U (q ) = 

27re 2 /(e^/g 2 a . + + k 2 ) is the Fourier transform of the 

screened impurity potential U(r) — (e 2 /47re)(e~ fcsr /r), 
where k s is the inverse screening length and e is the di- 
electric constant of the material. In the limit of small 
|q | <C k s , U(q a ) ~ 27re 2 /(efc s ) = Uq. Here, r and R are 
the position vector of electron and impurity, respectively. 
In this limit, we can use r 2 w ttIqIi 2 /NjIIq , where Ni is 
the 2D impurity number density. 

We follow Refs££ to calculate the collisional conduc- 
tivity. We include the correction to the unperturbed 
eigenstate \^f k (r)) due to the weak perturbative term 
AH = Hi + H2 + H3 + H3 . The first-order correction to 
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the Landau state is obtained by 

l^yrfc _ |, T ,± 



E 



(*t ikv (r)\AH\*t, tkv (r)) 
E x — E x i 



Following Refs.— £ and using the the perturbed Landau 
states |^^t (r))', we obtain the collisional conductivity 
as 



e 2 AWp 2 
h 27raFn 



± 7±1 



i?„ ,7, 



The exact expressions of if 1 and are given as 
If = {(2s T l)Dt-2sD 2 s +2s±l}/A 2 s , 

dky - 



M± = 



Jt = 



r a / l l 


r df{Ey 




'o 


l dE - 


E 


r a / l o 


r 
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E 



E=E±, 



E=E±,„ 



(30) 

(31) 
(32) 



sin (qxo)dk y . (33) 



The term Rf is appearing due to the first-order correc- 
tion to the Landau wave function and it is of the order of 
Vg . We neglect this term because of small contribution. 
The major contribution to the collisional conductivity is 
due to the first term proportional to Mf. The effect of 
the magnetic modulation mainly enters through the en- 
ergy correction due to the modulation in the total energy 
in the Fermi-Dirac distribution function. 

Similar to the diffusive conductivity, the collisional 
conductivity in the presence of the modulation will have 
two contributions, namely the SdH and the Weiss con- 



tributions: (f*L = atr 



SdH + ^Weiss- Jt is difficult to get the 
analytical expression of the Weiss contribution (c^, ciss ) 
comes from the energy correction in Fermi-Dirac distri- 
bution function. 

The numerical results of the change in the collisional 
conductivity A.a xx versus magnetic field for a = 2ao are 
plotted in Fig. (3). To compare the results of the mag- 
netic modulation case, we present A.a xx for the electric 
modulation case in the lower panel of Fig. (3). 

The oscillatory behavior with beating pattern appears 
in the changes in the conductivity due to the modulation 
at low magnetic field range where modulation strength is 
not much less than the energy scale of the Landau levels. 
The effect of the modulation diminishes with the increase 
of the strength of the perpendicular magnetic field. As 
the magnetic field increases, the SdH oscillation starts to 
dominate over the modulation induced Weiss oscillation. 

The collisional conductivities for spin-up and spin- 
down electrons oscillate with the same frequencies as the 
bandwidth | | . The frequencies of the oscillation of the 
bandwidths are the same as that of the diffusive conduc- 
tivities, namely, . Therefore, beating condition for the 
collisional conductivity is the same as given in Eq. (|2"Tj) 
for the diffusive conductivity. Using the beating condi- 
tion, we get the magnetic fields (in Tesla) corresponding 
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FIG. 3: Plots of the modulation induced change in the colli- 
sional conductivity versus magnetic field for a — 2an- 



to the beating positions, obtained from the Eq. (|2TJ) . are 
0.0506, 0.0574, 0.0662, 0.0782, 0.0957.... These are in 
good agreement with the exact numerical results shown 
in Fig. 3. The phase difference between the oscillations 
in the diffusive and collisonal conductivities is nearly tt. 

The modulation gives very small effect on the SdH part 
((7g dH ) of the collisional conductivity. It is also difficult to 
get an analytical expression of the SdH oscillations super- 
posed on the Weiss oscillations. Ignoring the modulation 
effect and using the analytic form of the DOS, the asymp- 
totic expression of the SdH oscillation has been studied 
in RefM. 

The beating pattern in the SdH oscillation is given by2£i 



E, 



1 + 



(To 4(w t) 2 

x cos(2tt fa/ B ) cos(2tt/ ' d /B 



(34) 



where cto = n e e 2 r/m* is the classical Drude conductivity, 
E F = [1 + E a /(2E F ) T (3/2)^E~/E7] and /„ = (/+ + 
D/2 and /«, = (/+- f-)/2. 
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C. The Hall conductivity 

The off-diagonal elements in the conductivity tensor 
are termed as the Hall conductivity which is given by^£ 



'yx 



ie 2 h 



£/e(i-/e)<£kim'H£) 



(35) 



To simplify the above equation, we shall use the following 
relation 



hO-fe) 



1 _ e P(E ( -E e ) 



k-k 



(36) 



The matrix elements of the velocity operators are zero 
except between the nearest Landau levels. Following 
Refs£&2!, we calculate the velocity matrix elements and 
the energy difference between two nearest Landau levels 
(see Appendix B). Finally, we get the expressions of the 
Hall conductivity as 



S 



1 



A S A S+1 

2 r a / l a 2 
J() 



D s+ i (D s ^/s ± V2k a lo 



[1 + cosqx 



-dk y . 



(37) 



In the large s limit, reduces to 
V B V\(p F TP a ) 3/2 



7TC, 
X COS 



sin {n/(p F TP J} 



{2n\(p F TPa) — — 



(Pf T P a 



;}■ (38) 



When T -> and E s < E F < E s+1 , Eqs. §7$ and 
give the asymptotic form of the Hall conductivity as 



^ [1 + 3 (T ± )2] 



where 



1 



1 (ka\ 2 

2\k F J . 



1 - 



h \ 2 

"-a 

kp 



(39) 



(40) 



is the total Hall conductivity of the 2DEG with Rashba 
SOI but without modulation and cr H = n e e/BQ is the 
classical Hall conductivity. Note that n e is the sum of 
the density of the spin-up and spin-down electrons. By 
taking the superposition of the analytical expression of 
the conductivity for the spin-up and spin-down electrons, 
we get the same beating condition as given in Eq. (|2"Tj) 
for the diffusive conductivity. The numerical results of 
the change in the Hall conductivity versus magnetic field 
is shown in figure 4. 

Similarly, we obtain the Hall conductivity for the elec- 



trical modulation case as a. 



-^[l + (3/2)(T± )% 
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FIG. 4: Plots of the change in the Hall conductivity due to 
modulation versus magnetic field Bo for a = 2qq. 



x10"' 




FIG. 5: Plots of the change of the Hall conductivity versus 
magnetic field So for different values of a. We have scaled up 
the electrical modulation case by 10. 



where 



V E 



X 



■sin{27r/(p F TP J} 



ne a V P F T Pa 
x sin|27rA(p F TPa) ~ tt/ 4 }- 
In figure 5, we compare the change in the Hall con- 



(41) 
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ductivity of the electric and magnetic modulation case. 
The beating condition remains the same in the electric 
modulation case also. The amplitude in the change in 
conductivity due to the magnetic modulation is nearly 
[p f /(2-k)] 2 times higher than the electrical modulation 
case. There is a tt phase difference between the Hall 
conductivities of the electrical and magnetic modulation 
cases. In the above analytical expression modulation ef- 
fect through the Fermi-Dirac distribution function has 
been ignored. 
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Appendix A 

We derive the DOS by taking the imaginary part of 
the self-energ y 39 ' 40 which is given as 



IV. CONCLUSION 



£-(£)■ 



(Al) 



We have studied magnetotransport properties of the 
2DEG in the presence of the Rashba SOI when the per- 
pendicular magnetic field is weakly modulated. 

The diffusive conductivity shows beating pattern due 
to the interference between the conductivities for spin- 
up and spin-down electrons. We calculate the asymp- 
totic expression of the Weiss conductivity which matches 
very well with the exact numerical results. The number 
of Weiss oscillations between any two successive nodes is 
fixed for given values of n e and a where as in the SdH it 
depends on the magnetic field. We have shown that the 
Rashba SOI strength can be determined by analyzing the 
beating pattern in the Weiss oscillation. The strength of 
the Rashba SOI can be determined by just counting the 
number of Weiss oscillations between any two successive 
beat nodes. There is a definite phase difference between 
the conductivities for magnetic and electric modulation 
cases. In magnetically modulated system, there is a dif- 
ference in amplitudes by a factor p a /(2ir 2 ) between the 
conductivity due to the spin-up and spin-down electrons 
where as in the electrically modulated system, the ampli- 
tudes are the same for spin-up and spin-down electrons. 

To observe the effect of the modulation, we plot the 
change in the collisional conductivity due to the mod- 
ulation at the low range of the magnetic field. It is 
found that the effect of the magnetic modulation is much 
higher than the electrical modulation. The major effect 
of the modulation comes through the energy correction in 
the total energy in the Fermi-Dirac distribution function. 
The beating condition in the Weiss contribution to the 
collisional conductivity is same as that of the diffusive 
conductivity. 

The modulation effect on the Hall conductivity is 
shown by plotting the change in the Hall conductivity 
due to the modulation. The beating pattern appears in 
the Hall conductivity and it increases with the increase 
of the Rashba strength. The oscillations are out of phase 
between the electric and magnetic modulation cases. The 
beating condition remains same as the diffusive conduc- 
tivity. The amplitude of the fluctuation in the presence 
of the magnetic modulation is found to be much higher 
in comparison to the electric modulation case. 



The DOS is the imaginary part of the self-energy: 
D(E) = Ira 



First we consider the lower branch. 



To find the summation, we use the residue theorem 
and neglecting (E a /frujo) 2 term, we obtain 'E~(E) — 

l gcot(7r S+ ) , where s + ~ ^{E - E~(E) + E a /2 + 

^E 2 + E a E}. We write E~(E) = A + iT/2, then the 
above equation can be re-written as 



A + ?T/2 = T-^cot 



iv) 



/ttT 

V huiQ 



2, 



sin it + i sinh v 

cosh v — cos u 
(A2) 



Here, u 



■{E - A + E a /2 + ^E 2 + E a E} and v 



(ill). 



TrT/(Hu ). The imaginary part is 5 - yj— j^- v _,^„ 
We are using the following standard result to simplify it 
further: 



sinhf 



cosh v — cos u 



= 1 



fc=i 



e kv cos(feu). 



(A3) 



We are considering the first term corresponding to k = 1 
only. Other terms are very small compared to the first 



one. We have ~ = (^j^j 



1 + Ze-^K^ cos(u) 



In 



the limit of 7rr 3> fiajQ, after first iteration, we have 
r/2 = ttTq/ (huo). We are putting it back in the ear- 
lier expression, we get 



r 

2 



1 



2 exp | 



cos {E + E a j2- 



( huQ ) i 



I Hujq 



El 



E n 



'«£)}]. (A4) 



Similarly, one can do it for the upper branch. Finally, the 
DOS for lower and upper branches can be put together 
as 



D ± {E) 



X cos 



2ith 2 
f 2;r 



foul® \ 2 



El 



E„E 



}jA5) 



In absence of the Rashba SOI and the Zccman term, in- 
cluding the spin degeneracy the above expression can be 



10 



reduced to the following standard resul 



39,40 



as 



where X 



■xq. Similarly, for the velocity component 



^^[ 1 + 2eXP {- 2 (^) I 



/2ttE 



cos 



V Hun 



(A6) 



Appendix B 

The velocity operators are given by 



dHo 

dp x 



Px 



a 

m* h 



(Bl) 



s+l,k 



, (0 I V V 



s , k 



y/A a A a 



2m* 



D S+1 (D S ^ ± V2k a l ) + V~s + pi.4) 



dH 
dp 



1 a 
— {Py + eB x) + —a x 



in 



Py+eB x 



m* 

n 



h 

Py+eBgx 



(B2) 



For £' : {s + 1, ±, k y }, the velocity matrix elements (see 
Refe^Sl) are 



-l,k. 



(r)) 



V, 



V A S A S 



D S D 



s^s+l 



(0f_x(X)|^|^(X)) 
m 

(<Pt(X)\^\4>f +1 (X)) Ti^D s+1 



VA7A 



s+l 



ftjUJn 



2m* 



+ Vs+T 



D s+1 (D s ^±V2k a l ) 
(B3) 



The energy difference between the two successive Landau 
levels is 



Es+i.k v - -fiwo[l - ~ff cos(gxo)], 



(B5) 



where = (F s+1 — F^)/(huo)- Substituting the above 
three equations (|B3|) . (|B4I) and (|B5[) in the Hall con- 
ductivity expression, we get Eq. (I37|) . A multiplication 
factor 2 need to be used for the contribution coming from 
£' : {s — 1, ±, fc y } states. 
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